Slow blow-up solutions for the H^1(R^3) critical focusing semi-linear wave equation in R^3

SLOW BLOW-UP SOLUTIONS FOR THE ff^R 3 ) CRITICAL
FOCUSING SEMI-LINEAR WAVE EQUATION IN M 3

J. KRIEGER, W. SCHLAG, AND D. TATARU

Abstract. Given u > i and <5 > arbitrary, we prove the existence of energy
solutions of

(0.1) d tt u-Au~u 5 =0

in R 3+1 that blow up exactly atr = i = 0ast-»0- . These solutions are
radial and of the form u = X(t) 2 W(X(t)r) + rj(r, t) inside the cone r < t, where
X(t) = t" x ~" , W(r) = (1 + r 2 /3)~i is the stationary solution of JOT}, and T}
is a radiation term with

/

J\r

(|Vf?(ai, t)\ 2 + \r) t (x, t)\ 2 + \ v (x, t)\ 6 ) dx^O, t->0

r<t]

Outside of the light-cone there is the energy bound

/ (|Vu(x,t)| 2 + \u t (x,t)\ 2 + |u(a:,i)| 6 ) dx < 5
J[r>t]

for all small t > 0. The regularity of u increases with v. As in our accompa-
nying paper on wave-maps [101 . the argument is based on a renormalization
method for the 'soliton profile' W(r).

1. Introduction

Since the seminal paper of Jorgens [5] much work has been devoted to the study
of well-posedness of the nonlinear wave equation

d tt u - Au + f(u) =

in R 3 ^ 1 and suitable nonlinearities f(u). Jorgens showed that for H 1 (R 3 ) subcriti-
cal defocusing nonlinearities f(u) = |u| p_1 u withp < 5 smooth data lead to smooth
solutions for all times. The critical defocusing casep = 5 was resolved by Struwe [H]
for radial data and Grillakis [S] for general data. These authors proved global well-
posedness and scattering results for energy solutions, see Shatah-Struwe [17] and
Sogge [15] . No corresponding results are known for the supercritical case p > 5.

In this paper we address the solvability of the nonlinear wave equation in M 3+1
with a focusing nonlinearity f(u) = — |u| p_1 m. In this case blow-up may occur.
Indeed, it was shown by Levine [llj via a convexity argument that data in (H 1 n
L 2 )xL 2 with negative energy lead to finite-time blow-up, see also Strauss [50] . Local
well-posedness in the optimal regularity class was considered by several authors, see
Sogge [18] for a detailed exposition of this work. Most relevant for us is the case
p = 5 where the equation is locally well-posed in the energy space H 1 x L 2 (R 3 ).
Moreover, if the solution cannot be continued beyond some finite time T* as an

The authors were partially supported by the National Science Foundation, J. K. by DMS-
0401177, W. S. by DMS-0617854, D. T. by DMS-0354539, and DMS-0301122.

1

2

J. KRIEGER, W. SCHLAG, AND D. TATARU

energy solution, then necessarily the Strichartz norm ||u||i8([ 0! T,)xR 3 ) — 00 (with
similar results in all dimensions).

The question of the blow-up rate was addressed by Merle-Zaag in the conformal
range p < 3, see jT3]— [TS] (their results extend to all dimensions). They showed
that if solutions to the Cauchy problem, with 1 < p < 3,

d tt u - Au ~ = 0, (it , ui) G Hl oc x Lf oc

blow up in finite time T*, then the following holds: for any a£l 3 the self-similar
change of variables

u(x, t) = (T* - ty^w a {y, s), y = ^ — - s = - log(T* - t)

leads to functions w a satisfying

sup \\w a (s)\\ H i(B) + \\d s w a (s)\\ L 2 {B) < K

s>- logT,+l, aeM 3

where B is the unit ball and a constant K that only depends on p, T* and the norm
of the initial data in Hl oc x Lf oc .

For the energy critical case p — 5, i.e.,

(1.1) d tt u-Au-u b =

there has been some recent activity, see [9] , [8] , [7] , which we now describe in more
detail. The Talenti-Aubin solutions

W(r) = (1 + r 2 /3)-^

are extremizers of the Sobolev imbedding if 1 (R 3 ) =— ► L 6 (M 3 ) and satisfy the Euler-
Lagrange equation —AW — W 5 = 0. In [5] the first two authors showed that there
exists a small co-dimension one manifold M. around W in a suitable topology so
that data on this manifold exhibit global existence and an asymptotic behavior of
bulk term plus radiation. The radiation term is also shown to scatter like a free
energy wave. It is conjectured, see [2], that this manifold has the property that
it separates a region of scattering from one of blow-up. As a first result in this
direction, Karageorgis-Strauss [7] showed that above the tangent space of M. at W
finite time blow-up occurs, albeit for the equation

d tt u-Au- \u\ 5 =

Note that the result of [9j equally well applies to the nonlinearity |u| 5 (in fact, the
solutions constructed in [9] are positive so that there is no distinction between u b
and I it 1 5 from the point of view of that paper).

Kenig-Merle [8] studied the behavior of solutions with data {uq,Ui) G H 1 x
L 2 (M 3 ) of energy S(uq,Ui) < £(W, 0) where the conserved energy is

S(u,u t )= ( [i(^ + |Vu| 2 )-^'
J R 3 L2 6 .

d.r

They found that in this regime there is a dichotomy between blow-up and global ex-
istence/scattering depending on whether ||Vuo||2 > ||VW||2 or ||V«o||2 < HVWlj 2-
Note that

W(x,X) := X^W(Xx)

SLOW BLOW-UP SOLUTIONS

3

is a stationary solution of (|1.1|) for all A > 0. Moreover, the energy is constant in
A (reflecting the energy criticality of the equation). Linearizing the wave equation
around W leads to the linearized operator

H = -A - 5W 4

The wave evolution of H has two types of instabilities: an exponential instability
arising from the negative spectrum of H (which has a unique negative eigenvalue)
as well as a "bound state" at zero energy: H(d\W\\=i) — where d\W decays
like r _1 and thus does not belong to L 2 (K. 3 ) — this is what one refers to as a zero
energy resonance. In this paper we construct blow-up solutions by 'projecting out'
the exponentially growing mode of the linearized equation.
More precisely, we seek radial, real- valued, blow-up solutions

u(x,t) = A(f)*W(A(t)x) +r)(x,t)

of where A(t) — > oo as t — > and with the local energy inside the light-cone
\x\ < t of T](x,t) going to zero as t — > 0. The local energy relative to the origin is
defined as

£ loc (v)= I (v? + \Vv\ 2 + \v\ 6 )dx

J[\x\<t]

The following theorem is the main result of this paper. The blow-up occurs at time
t = when solving backwards in time.

Theorem 1.1. Let v > ^ and S > 0. Then there exists an energy solution u of
(jl.ip which blows up precisely at r = t = and which has the following property:
in the cone \x\ = r < t and for small times t the solution has the form, with
\{t)=t- x - v ,

t) = \i(t)W(\(t)r)+ V (x,t)
where £\ oc (r](-,t)) — > as t — > and outside the cone u(x,t) satisfies

[ [\Vu(x, t)\ 2 + \u t (x, t)\ 2 + \u(x, t)\ 6 ] dx < 5

J[\x\>t]

for all sufficiently small t > 0. In particular, the energy of these blow-up solutions
can be chosen arbitrarily close to £{W, 0), i.e., the energy of the stationary solution.

The restriction v > 1/2 arises only due to technical reasons, and we hope to
eliminate it in subsequent work. If v > 1, then the solutions from Theorem 11.11
belong to L°°(R 3 ) for alH > and blow up at the rate

IM,f)||oo*H 1+ "V2

as t — > 0. The proof is based on a renormalization procedure analogous to the
one that the authors used for the construction of blow-up solutions for wave maps
in [10] . For our purposes this refers to the fact that we do not simply perturb around
X^(t)W(X(t)r) to obtain the linearized equation for -q, but rather first modify the
blow-up profile and then perturb around this " rcnormalizcd" profile. More precisely,
fix a large integer N. Then there exists a function u e satisfying

(1.2) u e e C*¥-({to > t > 0, |x| < £}), £i oc (u e ){t) < (tX(t))- 1 as t ^
so that the radiation term rj above has the form

rj(x,t) = u e (r,t) +e(r,t), 0<r<t

4

J. KRIEGER, W. SCHLAG, AND D. TATARU

where e decays at t = 0. In fact, r\ can be extended globally with the property that
£e^-(l 3 ), s t € ^-^-(M 3 ), £ioc(e)(*) <** as ^0

with spatial norms that are uniformly controlled as t — » 0.

As this paper was written concurrently with our wave-map paper |10] it is only
natural that there would be some similarities between this paper and [10] . In fact,
a secondary goal here is to show that the method used in both papers is flexible
and applies to quite distinct scenarios. The main differences between this paper
and [10] are as follows:

• The blow-up profile is not constant in L°° but rather grows at rate A 3 . The
renormalization procedure thus needs to be adapted to this case.

• In contrast to [10] , the linearized operator exhibits negative spectrum. This
produces exponential instability of the linearized wave flow.

• The linearized operator no longer exhibits a strongly singular potential in
the sense of [4]. Thus, a (Dirichlet) boundary condition is needed at R = 0.

We feel that the most important difference listed here is the exponential instabil-
ity. In fact, as in the asymptotic stability paper [9], we need to 'project out' this
exponential growth. Our blow-up rates are therefore expected to be non-generic.

2. The renormalization step

In this section we show how to construct an arbitrarily good approximate radial
solution to the wave equation (|1.1[) as a perturbation of a time-dependent ground
state profile

u = X^W(R), W(R) = {1 + R 2 /3rK R = rX(t)
with the polynomial timescale

A(t) = r 1 -", v > 0.

Theorem 2.1. Let k G N. There exists an approximate solution U2k-\ for
of the form

u 2k ^(r, t) = As ( t ) [W(R) + ^i? 2 (l + i?T 1/2 + O ( ^ 2(1 ( +jr r
so that the corresponding error has size

e2k~

1 yt 2 {txf k j

Here the O(-) terms are uniform in < r < t and < t < to where to is a fixed
small constant.

Remark 2.2. The u e in Q is

The analysis below shows that it has the stated regularity up to the light-cone.
Moreover, one checks that

A^

SLOW BLOW-UP SOLUTIONS

5

which is the claimed decay rate for the local kinetic energy of u e . The local potential
energy of u e decays like (iA)~ 3 .

Proof. We iteratively construct a sequence u k of better approximate solutions by
adding corrections v k ,

Uk = v k + Ufc-i

The error at step k is

2

e k = {-d 2 + d 2 + -d r )u k + u%
r

If u were an exact solution, then the difference

e = u — tifc-i

would solve the equation
2

(2.1) (-9 4 2 + ^ + -a r )e + 5^_ 1 £ + 10u|_ 1 e 2 + 10^_ 1 e 3 + 5ii fe _ie 4 +e 5 + e fc _i =

r .....

In a first approximation we linearize this equation around e = and substitute
Uk-i by uq. Then we obtain the linear approximate equation

(2.2) (^d 2 + d 2 + ^d r + 5u{\ e + e fc _! «

For r<fwe expect the time derivative to play a lesser role so we neglect it and
we are left with an elliptic equation with respect to the variable r,

(2.3) ^a 2 + ^ I . + 5u^e + e fe _i«0, r<£

For r w t we can approximate Up by zero and rewrite (|2.2p in the form

(2.4) ^-d 2 + a 2 + ^e + e fc _ 1 «0

Here the time and spatial derivatives have the same strength. However, we can
identify another principal variable, namely a — r/t and think of e as a function of
(i, a). Later, we reduce the above equation to a Sturm-Liouville problem in a which
becomes singular at a = 1.

The above heuristics lead us to a two step iterative construction of the w^'s. The
two steps successively improve the error in the two regions r <gC t, respectively r«t
To be precise, we define v k by

(2.5) (d 2 r + \ + 5u£j v 2k+ i + e° 2k =
respectively

(2.6) (-d? + d 2 + ^d r ^V2k+e° 2k _ 1 =0

both equations having zero Cauchy dat£0 at r = 0. Here at each stage the error
term e k is split into a principal part and a higher order term (to be made precise
below) ,

e k = el + el

lr The coefficients are singular at r = 0, therefore this has to be given a suitable interpretation

6

J. KRIEGER, W. SCHLAG, AND D. TATARU

The successive errors are then computed as

e 2k = e\k-\ + N 2k{v2k), e 2 fe+i = e\ k - d 2 v 2k+x + N 2k+ i(v 2k+ i)

where

(2.7) N 2k+1 (v) = 5(u 4 2k - -4) v + 10u 3 2k v 2 + I0u 2 2k v 3 + hu 2k v 4 + v 5
respectively

(2.8) N 2k (v) = f>u\ k _ r v + lOul^v 2 + 10u^_ 1 « 3 + hu 2k ^v 4 + v 5

To formalize this scheme we need to introduce suitable function spaces in the cone

C = {{r,t) : 0<r<t, < t < t Q }

for the successive corrections and errors. We first consider the a dependence. For
the corrections Vk we set

v — 1 1

and use

Definition 2.3. For i G N we let = if is is irrational, respectively j(i) = i if
v is rational.

a) For any positive integer k, we define Q to be the algebra of continuous func-
tions q : [0, 1] — ► K with the following properties:

(i) q is analytic in [0, 1) with an even expansion at 0.

(ii) Near a = 1 we have an absolutely convergent expansion of the form

oo

q(a) = q (a) + £(1 - af^~ 2 ^ £ ««(o)(Iog(l - a)y

i=l ]=0

with analytic coefficients q , qij.

b) Q m is the algebra which is defined similarly, with the additional requirement
that

3ij(l) = if i = 4fc + l > 4m+l.

We remark that the exponents of 1 — a in the above series are all positive because
of Po > —\- For the errors we introduce

Definition 2.4. Q' is the space of continuous functions q : [0, 1) — > M. with the
following properties:

(i) q is analytic in [0, 1) with an even expansion at

(ii) Near a = 1 we have a convergent expansion of the form

oo

q(a) = q (a) + ]T(1 - af^- 2 ^ 1 ]T gy (a)(log(l - a)Y

i=l 3=0

with analytic coefficients qo, q^ .

b) Q' m is the space which is defined similarly, with the additional requirement
that

q^ (1) = if i = 4k + 1 > 4m + 1.

SLOW BLOW-UP SOLUTIONS

7

By construction, Qk C Q' k . The families Q! and Q' k are obtained by applying
a~ 1 d a to the algebras Q and Qk, respectively.

We remark that the number of logarithms in these definitions in the case when
v is rational is far from optimal, but we have chosen this form since it simplifies
the presentation. Next we define the class of functions of R:

Definition 2.5. S m (R k ) is the class of analytic functions v : [0,oo) — ► K with the
following properties:

(i) v vanishes of order m and R~ m v has an even Taylor expansion at R = 0.

(ii) v has a convergent expansion near R — oo,

k-2i

i=0

The importance of even expansions in R lies with the fact that only those corre-
spond to smooth functions in R 3 . For the same reason, we will work with even m.
We also introduce another auxiliary variable,

^ 6 =w

Since we seek solutions inside the cone we can restrict b to a small interval [0, bo]-
We combine these three components in order to obtain the full function class which
we need:

Definition 2.6. a) S m (R k , Q n ) is the class of analytic functions v : [0, oo) x [0, 1] x

[0, bo] -> K so that

(i) v is analytic as a function of R,b,

v : [0, oo) x [0, bo] -> Qn

(ii) v vanishes of order m and R~ m v has an even Taylor expansion at R = 0.
(Hi) v has a convergent expansion at R — oo,

v(R,;b) = Y,Ci(;b)R k ~ 2i

i=0

where the coefficients Ci : [0, bo] — > Q m are analytic with respect to b.

b) lS m (R k , Q n ) is the class of analytic functions w on the cone Co which can be
represented as

w(r,t)=v{R,a,b), veS m (R k ,Q n )

We note that the representation of functions on the cone as in part (b) is in
general not unique since R, a, b are dependent variables. Later we shall exploit this
fact and switch from one representation to another as needed. We shall prove by
induction that the successive corrections Vk and the corresponding error terms
can be chosen with the following properties: For each k > 1,

8

J. KRIEGER, W. SCHLAG, AND D. TATARU

(2.10) V 2k -1 G T^fe K 2 (^, Qk-l)

(2.H) t 2 e 2k ^ 1 e 77 ^lS Q (R,Q' k _ 1 )

(txy
(txy

(2-12) ^G^y^IS 2 ^ 3 ^)

(2.13) t 2 e 2k G j^k-[lS (R-\Q k ) + blS (R,Q' k )]

with (|2.13[) also valid for k = 0. We remark that the order of vanishing at R =
can be successively improved with k, but this does not appear to be important.

Step 0: The analysis at k =

With our notations, one checks that

(2.14) t 2 e = -t 2 dtt[X^W{X(t)r)] G X^IS°(R- 1 )

as claimed. Now assume we know the above relations hold up to k — 1 with k > 1,
and we show how to construct v 2k -i, respectively v 2kl so that they hold for the
index k.

Step 1: Begin with e 2k - 2 satisfying (|2.13[) or (|2.14[) and choose v 2k _\ so that
(f2~T0l) holds.

If fc = 1, then define e[J := eo- If k > 1, we use (|2.13|) to write

e2fc-2 = e 2fe _2 + e 2fc-2

where

* e K-2 G 7^Z2 IS°(iT \ Q fc _x), i 2 e 2fe _ 2 G ^fe IS°(iZ, Q' fe _ x )

In the first term we can set b — and eliminate the b dependence, as all the b
dependent part can be included in the second term.

We note that the term e 2fc _ 2 can be included in e 2 fc-i, cf. (|2. 1 1[) . We define
v 2k -i as in (|2.5[) neglecting the a dependence of e 2fc _ 2 . In other words, a is treated
as a parameter. Changing variables to R in (|2.5p we need to solve the equation

(tX) 2 Lv 2k ^ = t 2 e° 2k _ 2 G fc _ 2 IS ^" 1 , Q fc _x)
where the operator L is given by

L = -d 2 R - ld R - 5W 4 (R)
Then (|2. 10|) is a consequence of the following ODE lemma.
Lemma 2.7. The solution v to the equation

Lv = / G S°{R- 1 ), v(0) = v'(0) =

has the regularity

v G S 2 (R)

SLOW BLOW-UP SOLUTIONS

9

Proof. Since / is analytic at with a constant leading term, one can easily write
down an even Taylor series for v at with a quadratic leading term.

It remains to determine the asymptotic behavior of v at infinity. For this it is
convenient to remove the first order derivative in L (to achieve constancy of the
Wronskian). Thus, we seek a solution of

LR v = Rf, L = d% + 5W 4 = d 2 R

(l + i? 2 /3) 2

We use this fundamental system of solutions for L:

4>{R) = R(l - i? 2 /3)(l + i? 2 /3)"3

9(R) = (1 + i? 2 /3)^(l - 2R 2 + i? 4 /9)

Clearly, Ld\W = and we set <f> = Rd\W\\ = \ up to a constant. The function
9 is then determined from the Wronskian constancy condition W(6, cf>) = 1. This
allows us to obtain an integral representation for v using the variation of parameters
formula, which gives

[■R rR.

v = —R~ 1 9(R) / cj)(R')R' f(R!) dR' + RT^R) / 9(R')R'f(R') dR'
Jo Jo

The right-hand side grows like R, as claimed. □

(2.15) Vl= v^R), VeS 2 (R)

As a special case of the above computation we note the representation for v\ ,

(tXf

This justifies the choice of the second term in the expansion for it2fc-i in Theo-
rem [2TTJ

Step 2: Show that if V2k-i *s chosen as above then (|2.1ip holds.

Thinking of i^fc-i as a function of t, R and a we can write e2k-i in the form

e 2 fe-i = N 2 k-i(v 2 k-i) + E t v 2 k-\ + E a v 2k -i

Here N 2 k-i(v2k-i) accounts for the contribution from the nonlinearity and is given
by (|2.7[) . E t V2k-\ contains the terms in

(2.16) -d tt v 2 k-i(t,R,a)

where no derivative applies to the variable a, while E a V2k-i contains those terms

(d tt - rr - ^dr^j U 2 fc-1 (*, R, a)

r

where at least one derivative applies to the variable a (recall that in Step 1 the
parameter a was frozen). We begin with the terms in Nzk—i- We first note that,
by summing the Vj over 1 < j < 2k — 2,

(2.17) u2 k -2-u £ -^lS 2 (R,Qk-i)

(txy

The first term in iVafc— i(u2fe— l) contributes

t 2 ( u 2k-2 ~ U o) V 2k-l = < 2 [(u 2 fe-2 - Mo) 4 + 4(u 2 fe-2 - «o) 3 Wo

(2.18) +6(u 2 fc- 2 - u ) 2 ul + 4(u 2/£ - 2 - u )ul}v 2 k-i

10 J. KRIEGER, W. SCHLAG, AND D. TATARU

Using (|2.17|) we compute

A^

C a 6 IS 2 (iT 2 , Q fe -i)77^ IS Z (-R, Q fe _i)
A*

c (^F Is2(irl ' Sfc -

as well as

t 2 {u 2k -2 - uo)u 3 v 2 k-i G t 2 ^ IS 2 (i?, Q fc _ 1 )Ai5°(i?- 3 )— ^IS 2 ^, Q fe -i)

The other two terms in (|2.18p are similar. Next, compute

and

A

with similar statements for W2fe-2 u 2fe-i an d M 2fc-2f2fe-i- Summing up we obtain

(|^ IS2(rl ' St - l)C (|

This concludes the analysis of A^fc-i^fc-i)- We continue with the terms in
E t V2k-i, where we can neglect the a dependence. Therefore, it suffices to com-
pute

^((^ IS2(i?) ) C (^ IS2(i?)
Finally, we consider the terms in E a V2k-i- With

(U) 2fe

we have

A 5 \ N A*

(OP J «-a(i?,«) + OT

-2Ra~ 1 w Ra - 2aT 1 w a (i?,a) + (a 2 - l)w aa (R, a) + 2aw a (R,a)

V2k-i(r,t) = , w(R,a), w G 5 (i?, Qk-i)

t 2 E a v 2k -i = -2tdt ( y-^2T aw a (i?, a) + — ^ [2(i/ + l)aJ?«; afl ( J R, a)

SLOW BLOW-UP SOLUTIONS

11

Since Qk-i are even in a we conclude that

(1 - a 2 )d aa , ad a , a~ 1 d a : Qk-i -> Q^-i

and therefore

\|

i^ 2fc _ lG _^IS 2 (i*, Ci_x)

This concludes the proof of (|2.11|) . We remark that for the special case of k = 1,
i.e., with «i as in (|2.15p . these arguments yield

(2.19) t 2 ei G j^IS 2 (R)

Step 3: Define v 2 k so that (|2.12[) holds.

We begin the analysis with e2k-i replaced by its main asymptotic component e°
around R = oo. This has the form

which we rewrite as

(2-20) t 2 4 k _, = 7^p9(a), q e QLi

A^

t 2 e»-i = (a)2fc -i "g( a )
We remark that (|2.19[) implies that

* 2 e~?(a) = ^a.

Consider the equation (|2.6p with eS^^ on the right-hand side,

2,
r

We look for a solution fi2fc which has the form

A^

t 2 ( -d 2 + d 2 + -8 r i 2k = -#e\ h _ x

V2k

Wr~' W2k{a)

Thus,

t 2 [-d 2 + d 2 + h^j U^W 2k (a)^ = —^aq(a)
Conjugating out the power of t we get

t 2 ( - (ft + ( 2fc -l>-^ 2 + 9 2 + Hft) ^ (o) = ag(a )

which we rewrite as an equation in the a variable,

(2.21) L (2fc-!)»-| W M = aq(a), q G
where the one-parameter family of operators Lp is defined by

(2.22) Lp = (1 - a 2 )d aa + 2(a- 1 + a/3 - a)9 a - (3 2 + f3

We claim that solving this equation with zero Cauchy data at a = yields a solution
which satisfies

(2.23) W 2 k{a) = a 3 qi (a), q x G Q fe

12

J. KRIEGER, W. SCHLAG, AND D. TATARU

This gives

v * k = (t\)2k-i a 91 (°) = ( tA )2fc+2 J? 91 (a)

which is not entirely suitable as the next correction since it has an odd expansion
at R — instead of an even one. To remedy this we simply set

- R3 (irW 9l(a) e (^ 54(i?3 ' Sfe)

Clearly, this will conclude Step 3. To prove the claim (|2.23p we need the following

Lemma 2.8. Let f G aQ' k l , k > 1. Then there is a unique solution w £ a 3 Qk to
the equation

(2.24) L (2k _s )v _iw = f, w(0) = 0, d a w(0) =
Proof. Denote

(2.25) /? = (2fc - - 1 = (4fc - 3)/3 + 2(fc - 1) > ~.
We write

= a- 2 d a (a 2 d a ) - a 2 9 QQ + 2(/3 - l)a9 - /3 2 + f3

To study the behavior of the solutions at we match the coefficients in Lp w = f
with

oo oo
3=1 3=2

yields the system, with j > 1,

(2j + l)(2j + 2)^+1 + [-(2j - l)(2j - 2) + 2(/3 - l)(2j - 1) + (/? - /3 2 )] Wj =

where we take i«i = 0. The coefficient of Wj+\ is always nonzero; this allows us
to successively compute the coefficients uij. The convergence of the series for w
follows from the convergence of the series for /.

It remains to study the solution w near a — I. The behavior of Lp at 1 is well
approximated by

L} 3 = 2(1 - a)d aa + 2(3d a = 2(1 - a)" +1 o [(l - a)~ fi d a ]

which annihilates the functions 1 and (1 — a)^ +1 . Therefore, we seek a fundamental
system for Lpy — of the form

oo oo

(2.26) &(o) = l + ^/i,(l-a) f , 2 (a) = (1 - af +1 [l + ^ fi e (1 - a) 1
This leads to the conditions, with no = (Xq = 1,

t-i

(2.27) 2nt +1 (t +l)(t-p) = hW -2(3 + 3) + (3 2 -P}+ 2 £ j H

3=0

(2.28) 2fi e+1 (e + (3 + 2)(£+l) = fi e [([3 + I + £)(£ - [3 + 4) + [3 2 - [3]

i-\

+ 2^ + 1+^

SLOW BLOW-UP SOLUTIONS

13

Clearly, (|2.28|) always has a solution whereas (|2.27p requires (3 £ Zj; in the latter
case, the series in (|2.26p define entire functions. If, on the other hand, (3 £ Zq ,
then cf>i is modified to

oo

(2.29) <t>i{a) = 1 + ^^(1 -a) £ + Cl 2 (a)log(l -a)

e=i

with some unique choice of c\.

Modulo a linear combination of 4>\ , 4>2 it suffices to find one solution to the
inhomogeneous equation Lpw = f near a = 1. We begin with the case when / has
the form

/( ) = (i- a y>XA(i-

with 7 > 0. Then we seek a solution w of the form

a)

w{a)

k=0

This leads to the system
(2.30)

2^o(7+l) (7-/3) =/o

2w e+1 ( 1 + i + 2)( 7 -13 + 1+1) =f t + «j € [( 7 + 1 + 1)( 7 + 1 + 4 - 2/3) + /3 2 - (3\

2e£S(7+j + i:

<7=0V i T J T

which is solvable unless (3 — 7 is a nonnegative integer, in which case the solution
w has the modified form

00

w(a) = (1 - a) 7+1 X>fc(l - + c0 2 (a) log(l - a)

Indeed, assume that j3— 7 = p > with integer p. Then (|2.30[) is in general violated
for 1+1 = p regardless of the choice of w p ; thus, c above has to be chosen so that
the coefficients of (1 — a) 7+p = (1 — a)' 3 in Lpw — f match. This can be done since

<h(a) log(l - a) = [(1 - af +1 + 0((1 - af+ 2 )} log(l - a)

and

L/s[0a(o) log(l - a)] = 2(/3 + 1)(1 - a)' 3 + 0((1 - a)' 3 "*' 1 )

with (3 + 1 > i and thus nonzero. Note that w p remains undetermined - this
amounts to the freedom of adding a constant multiple of 02(a) to w(a).
Similarly, if / has the form

3 j 00

/(a) = J2 (M 1 - a)) m fm(a) = ]T (log(l - a)) m (l - a)~< £ f km (l - a) k

m=0 m=0 fe=0

then we seek a solution u> of the form

3 3 00

w(a) = J2 ( lo g(! ~ a)) m Wm(a) = £ (log(l - a)) m (l - a) 7+1 £ «, fem (l - a) k

m=0 m=0 fc=0

Identifying the coefficients of the powers of log(l — a) we obtain the system
L Wj = fj, Lpwj-x = fj-i - jQiwj,

14

J. KRIEGER, W. SCHLAG, AND D. TATARU

respectively

LpWj^k = fj_i - jQiWj-k+i + j(j - l)(9 2 w-,- fc+ 2, k > 2

where

l + o 2(1 + a) 2a(3 „ 1 + a
1 — a a 1— a 1— a

This system is solved iteratively as in the first case provided that (3 — 7 is not a
nonnegative integer. Otherwise, the solution w has the modified form

j 00
w(a) = £ (log(l - a)) m [c m 2 (a) log(l - a) + (1 - a) 7+1 £ w km (l - a) k ]

m=0 k=0

Consider now / 6 Q'u-i- If v i s irrational then according to Definition 12.41 we
can represent it in the form

f(a) = f (a) + £(1 - aT^-^^Ma)
i=l

with fi analytic. Hence the exponent 7 above takes the values

~i - 1

7, = Ufa + 1) - 2

- 1

On the other hand, we have

P = (2k - |)„ - 1 = (4k - 3)(A) + 1) - 2k + 1 = 7 4fe-3

Then 7^-/3 can only be an integer if i = Ak — 3. However, in this case there is
no logarithmic term due to the additional condition Q4fc_3(l) = which has the
effect of replacing 74^-3 by 74/C-3 + 1. The conclusion of the lemma follows in the
irrational case.

On the other hand, if v is rational then / has the representation

00 i

f(a) = f (a) + £(1 - )*CA.+i)-*[*ii]-i £ / tf (o)(log(l - a))'

i=l j=0

with fij analytic. In order for — "fe to be a nonnegative integer we need to have
I < Ak — 3. Again the condition 54^-3(1) = guarantees that there is no extra
logarithmic term arrising from the i = Ak — 3 component of /. On the other hand,
if I < Ak — 3 then we can contribute an extra logarithm for a total of I + 1 < Ak — 3
logarithms to the i = Ak — 3 term of w. The conclusion of the lemma again follows.

It is also worth discussing the special case k = 1. This will also serve to explain
how the algebra Qk arises in the iteration. If k = 1, then (|2.21j) reduces to the
equation, with the usual /3 = (v — l)/2,

L/3o W 2 (a) = a

due to t 2 fi(a) = \i(t\)~ 1 a. As discussed above,
W 2 (a) = g (a) + 9l (a)(l - af« +1 if (3 & Z+

W 2 (a) = h (a) + hi(a)(l - a) r3o+1 + h 2 (a)(l - af 0+1 log(l - a) if O G
Thus, we see that in all cases W 2 G Q\ for j = 0, 1 and a near 1. □

SLOW BLOW-UP SOLUTIONS

15

Step 4: With v 2k as above show that e 2k is as claimed.
Modifying (|2.20[) to insure an even expansion at R = we set

2fc 1 (t\) 2k (R2 + l)i k 1

Then we can write e2/c in the form

k-l) + ^ (4k-! - (-9t + + 2 -d r ^j V 2k \ + f N 2k (v 2k )

t 2 e 2k = t 2 (e 2 fc_i - e

where A^/c is defined by

We begin with the first term in e 2k , which has the form

t 2 (e 2fe _a - 4 k _ t ) G IS°(7r\ e' k _i)

We claim that

(2.31) IS°(i?- 1 ,Q' fc _ 1 )c IS°(i?- 1 ,Q fe _ 1 ) + ^ I IS 2 (i?,Q^ 1 )

For iu G IS^iZ -1 , Qfc_ x ) wc write

tu = (1 — a 2 )w + R 2 w

{txy

Then

as desired.

Next, we consider the expression

/ = * 2 (^ fc -i-(-d 2 + d 2 + V) « 2fe

By construction this would vanish if the factor R(l + i? 2 ) _ 5 were dropped in both
e 2fc-i an d w 2fc- Hence / contains only the components of the second term for which
at least one derivative falls on the factor R(l + R 2 )~^:

2„\ R „9„_„ R

f = V 2k t Z ( -d 2 + 31 + -d r ) - ^ - 2t 2 d t V 2k d t

(l+i? 2 )3 (l+i? 2 )2

2t-d r v 2k d r - — T +v 2k — o r -

(l + i? 2 )2 r (l + i?2) 5

where t>2fc = v 2 kR{^ + R 2 )~^ ■ Then a direct computation gives

as needed.

Finally we consider the nonlinear terms in N 2k {v 2k ). Again the a, 6 dependence
is uninteresting since Q k is an algebra. We start with the last term in (|2.8p . By
construction,

A3 _ .

16

J. KRIEGER, W. SCHLAG, AND D. TATARU

Thus,
Using that

t 2 v 5 2k g IS 5 (i? 5 , Q k ) c {a)2k (u)8fc _ 2 (l + i? 2 ) 3 IS 1 ^" 1 , C*)

(U) 8fc - 2

we conclude that

(1 + R 2 f _ ^ + + 3a 4 6 4fc-3 + a 6 6 4fe-4

t 2 v* k €-^ IS 1 (R-\Q k )

which is an admissible contribution to (|2 . 13|) . Finally, we check the first term in
dUg]) . From

A

(txy

and the form of V2k,

X^ „o,„ A 4 A*

w 2 fe-i -wo£ 77777 IS 2 (i?, Qfe)

^ II (/,') + -—IS^Q^J ' 2fc+2 IS 2 (R\Q k )

C ^A^iT 1 ) + X^a 2 IS°(R-\ Qk)) i IS 2 (i? 3 , G fc )

c (a) 2 is°(i?- 4 , g fc ) * 2 2 fc+2 IS"! A )

The other three terms in (|2.8p can be checked similarly. This concludes the proof
of Theorem O □

3. The linearized problem

We seek a radial solution of fTTTj) of the form it = u 2 fc-i + e with U2k-i as in
Theorem 12. II Our ansatz leads to

(3.1) d tt e -Ae- bX 2 (t)W A {X{t)x)e = N 2k -i(e) + e 2k -i

where N 2 k-i(£) is as in (|2.7p but with u 2 k-2 replaced with u^k-i- Set e(t,x) =
v(r(t), X(t)x) and y = \(t)x, r = r(t). Then, with A = ^,

fte(t,r) = r'(t)(w r + AA" 1 ^)

and

(3.2) tt e(t, r) = r"(t)(9 r + X\~ l yd v )v + r'(t) 2 (d T + XX~ 1 yd y ) 2 v
Set

/•to 2 !

r(t) = / A(s rfs + -to v = -t~ v
Jt v v

so that t' (t) = A(t), and r"(i) = A(r)A(r) (we are writing A(r) instead of A(t(r))).
Then (|3.ip can be rewritten as

[{d T + XX^ydyfv + XX- 1 (d T + XX~ 1 yd y )v ~ Av - hW A v]{T : y)

= A- 2 (r)[7V 2fc _ 1 ( e ) + e 2fe _ 1 ](<(r),A- 1 y)

SLOW BLOW-UP SOLUTIONS

17

We remark that the linear Schrodinger operator H = —A — hW A on L 2 (R 3 ) has
at least one negative eigenvalue as well as a zero energy eigenvalue and resonance.
The negative spectrum renders the linear evolution in (|3 . 3[) exponentially unstable.
To address this further, we switch to the radial variable R = Xr from the previous
section. Thus, (13.31) is the same as

[(d T + XX~ X Rd R ) 2 v + XX- 1 (d T + XX-'Rd^v - v RR - \v R - 5W 4 v}{t, R)

= A- 2 (r)[^ 2fe _i( £ ) + e 2 fc-i](t(r), X~ 1 R)
or, with the new dependent variable e(r, R) := Rv(r, R),

(d T + XX~ 1 Rd R ) 2 e - AA _1 (d T + XX~ 1 Rd R )i + Ce
(,) ' 4 ' = X- 2 R{N 2k -i(R- 1 e) + e 2k - 1 ]

where

C = -dun - 5W i (R) on L 2 (0,oo)
with a Dirichlet boundary condition at R = 0. Let

0(t) := AA- X (r) = V := d T + P(r)Rd R

TV

and rewrite (13.41) in the form

(3.5) V 2 e - (3{t)V£ + Ce = f

To solve this equation we need precise spectral information on the operator C.

4. The spectral and scattering theory of the linearized operator
Definition 4.1. Let ^

C:= ~ dRR ~0TR 2 JW 2

be the half-line operator on L 2 (0,oo) with a Dirichlet condition at R = 0. It is
self-adjoint on the domain

Dom(£) = {/eL 2 ((0,co)) : /, /' e AC([0, R}) V R, /(0) = 0, /" 6 L 2 ((0,oo))}

Note that C<fi = where

cj)(R) := 2Rdx ^ \*W(XR) = R(l - R 2 /3)(l + R 2 /3)-^

This means that C has a resonance at zero energy. Since <j) has a single positive zero,
it follows from oscillation theory, see [3J, that there is an unique simple negative
eigenvalue which we denote by Thus, there is cj>d € L 2 (0,oo) n C°°([0, oo)),
decaying exponentially, and with (j>d(R) > for R > but (^d(O) = so that
£-4>& — £,d4>d- We also assume that H^dlla = 1- Clearly, C has no other eigenvalues
or resonances.

Lemma 4.2. The spectrum of C equals

spec(£) = {£ d }U[0,oo).

The positive spectrum is purely absolutely continuous, and ^ is a simple eigenvalue
with eigenfunction that we denote by 4>d- Moreover, C has a resonance at zero.
In fact, C(/)o = with (f>a(R) = R(l - R 2 /3)(l + R 2 /3)'i. Finally, C is in the
limit-point case at infinity.

18 J. KRIEGER, W. SCHLAG, AND D. TATARU

As usual, see Marchenko [12] or Section 2 of [4], we introduce the standard
fundamental system of solutions (f>(R, z) and 9(R, z) for all z G C of Cy = zy with
the boundary conditions

0(0, z) = 9'(0, z) = 0, <f>'(0, z) = 9(0, z) = 1

so that in particular

W(6(;Z),<t>(;Z)) = l

These functions are entire in z. Note that 4>(R, 0) = 4>o(R) from above. Further-
more,

9 (R) := 9(R,0) = (1 - 2R 2 + i? 4 /9)(l + i? 2 /3)"i
The Weyl-Titchmarsh function m(z) is uniquely defined by

(4.1) V+0) z ) : = z ) + ™>(z)4>(;z) £ L 2 (0, oo) V Imz >

The solution ip + is referred to as the Weyl-Titchmarsh solution. Then one has the

following, see H):

Proposition 4.3. The function m can be analytically continued to C\spec(£) and
it is a Herglotz function. For each R>0, ip+(R, z) and i/j + (R,z) are analytic on
C\ spec(£). The spectral measure of C is

dp = S Sd + p(0 = - Im m(C + t0)

7T

in the following sense: the distorted Fourier transform defined as

^ : / — ► /

poo pb

/(&)=/ MR)f(R)dR, /(0=lim / <t>(R,0f(R)dR, £>0
is a unitary operator from L 2 (K. + ) £o L 2 ({^} U R + , p) and its inverse is given by

t- 1 : f — ► /(i?) = ffa)<t> d (R) + lim /" W 0/(0 p(£) d£

Here lim refers to the L 2 (R + ,p), respectively the L 2 (R + ), limit.

In the sequel we view the Fourier transform as a vector-valued map

7(60

f

1 /(•)

The Weyl-Titchmarsh solutions are scalar multiples of the Jost solutions f+(R, z)
which are determined via the condition that

(4.2) Cf + (;z) = zf+(;z), f+(R,z)~e i ^ R BB R^>00

where Imz > 0, Im^/z > 0. They are solutions to the integral equation, with
V(R) = -5(l + R 2 /3)-\

f + (R, z) = + ^ Sin ^ {R '- R)) V{R')f + {R!, z) dR>

Jr V z

The functions ip+(R,z) have well-defined limits as z — > £ + iO when £ > 0. In

i

particular, tp + (R,£) ~ c (£)e l ^ R as R — » oo. The constant c (£) is determined
from the Wronskian condition W(?/'+(-, £)> 0) = 1- Once we have determined
c o(£) we find m(£ + «0) from the Wronskian relation

m(£ + »O) = W(0(-,0^+(-.£ + iO))

SLOW BLOW-UP SOLUTIONS

19

We now give an asymptotic expansion of our fundamental system for small z.

Proposition 4.4. For any z€C the fundamental system <j)(R,z), 6{R,z) admits
absolutely convergent asymptotic expansions

oo

<f>(R,z) = MR) + R- 1 J2( R2z ) j( t>AR 2 )

oo

0(R, z) = O (R) + J2( R2z ) j0 j( R2 )

i=i

where the functions cfrj, 9j are holomorphic in O = {ti 6 C : Reu > —5} and
satisfy the bounds

&
(~)!

(4-3) |0 3 -(u)| < —- — \u\(u)-*,

(4.4) <__( u ) 5 , ue n

Furthermore,
(4.5) ct>x{u)

-\u{l + o(l)) as

(l + o(l)) as moo

Proof. We begin with <fi. We formally write

<f>(R,z) = R- l Y,^fA R l MR) = R4>o(R)

3=0

This becomes rigorous once we verify the convergence of the series in any reasonable
sense. The functions fj should solve

CiR- 1 ^) = R^fj-i, fj(Q) = /j(0) =
where we have set /_i = 0. The forward fundamental solution for C is

h(r,r') = (MR)W) - MR')HR))1[r>r>]

Hence we have the iterative relation

f R R r i
= / 157 MR)W) ~ MR'Wo(R) fj-i(R')dR', f (R) = RMR)
Jo ft 1 J

Using the expressions for tfio, 9q we rewrite this as

fj(R) =

" R R 2 {\ - i? 2 /3)(l - 2R' 2 + i?' 4 /9) - RR'(1 - i?' 2 /3)(l - 2R 2 + i? 4 /9) . , D „ , D ,

5 -. fi—\{R ) dR

R'{l + R 2 / , Z)i{l + R' 2 /Z)i J

It is clear from this integral that each fj(R) is an analytic function of R 2 provided
Rei? 2 > —1. Moreover, fj(R 2 ) vanishes like i? 2 ^ +1 ^ around R = 0, and grows
like i? 2j+1 as R — > 00. The bounds in (14. 3p are a quantitative version thereof, and
proved by induction.

For 6(R, z) we make the ansatz

00

6(R, z) = J2 z J 9j(R), 9o(R) = o (R)
3=0

20

J. KRIEGER, W. SCHLAG, AND D. TATARU

where the functions gj solve

£9j=9j-i, .9-1=0

Thus, the iterative relation is

MR)0o(R') - MR'Wo(R) gj-i(R') dfi, g (R) = (R)

The analyticity is the same as for the fj(R), and gj(R) vanishes like R 2j around
R = 0, whereas the growth is R 2j+1 as R — ■> oo, as claimed.

Finally, the leading order of 4>\ (u) is found by solving for the coefficients c\ , C2

in

(-d RR -5(l + R 2 /3)- 2 ) Cl R 3 (l + o(l)) = R(l + o(l)) as R -►
(-<9fli?,-5(l + i? 2 /3)- 2 )c 2 i? 2 (l + o(f)) =-V3(l + o(l)) as i? ^ oo,
respectively. Thus, c\ = — g, c 2 = as claimed. □
Next, we turn to the asymptotic expansion of the Jost solutions.

Proposition 4.5. For any £ > 0, f/ie Josi solution /+(■,£) as m (|4.2|) is of the

form

f + (R,Z) = e iR ^o-(R^,R), R 2 t>l
where a admits the asymptotic series approximation

oo

a(q,R)^J2^^t( R )

3=0

in the sense that for all integers jo > 0, and all indices a, f3, we have

30

(4.6) sup (R) 2 (Rd R ) a (qd q f a(q, R) - T q- j 4>f(R)

R>0 1

3=0

< C,

c*Mo q

-i'o-l

for all q > 1 . Here

4'+ = l, 4t(R)

as R — ► oo
as fi^O

zcii?" 2 + iO(R- 4 )
ic 2 R + iO{R 2 )

with some real constants c\ , c 2 . More generally, ipt (R) o,re smooth symbols of order
—2 for j > 1, i.e., /or all k >

sup(R) 2 \((R)d R ) k 4UR)\<oo

R>0

Finally, ipf(R) = 0{R : >) as 0.
Proof. With the notation

O -( 9)J R) = / + (i?,0e- i ^
we need to solve the conjugated equation

5

(4.7)

-d RR - 2i^d R -

(l + i? 2 /3) 2

We look for a formal power series solving this equation

a(R^,R) =

(4.8)

J=0

SLOW BLOW-UP SOLUTIONS

21

This yields a recurrence relation for the fj's,
with /o = 1, which is solved by

tm - - 2 d R f^{R) - - 2 J r jYTww fj - l{R,) dR '

Then fi(R) is smooth for all R > with /i(-R) = -5^ + iO(R- 5 ) as i? -> 00 and
= ic + iO(R) around # = 0. More generally, = iW(R- j ~ 2 ) as i? -> 00

and fj(R) — VO{X) as i? — > 0. Differentiating these functions leads to symbol-type
behavior as R — » 00. Defining ip~j~(R) :— R 3 fj(R) yields the desired bounds.

To finish the proof, we first construct an approximate sum, i.e., a function
cr ap (q, R) with the property that for each jo > we have

(4.9) \(Rd R ) a (qd q f[a ap (q,R) - f>-^t(i?)] | < c a)/3j0 (i?)- 2 ^' 0-1

i=0

The construction of (J a p{q, R) is standard in symbol calculus; we set

00

a ap (q,R) :=5^g-^t(Ji)x(?«i)
i=o

where #j — > sufficiently fast and x is a cut-off function which vanishes around zero
and is equal to one for large arguments. The bound (|4.9[) implies that cr ap (i?^5, i?)
is a good approximate solution for (|4.7|1 at infinity, namely the error

e(i?^,i?) = (-fcut - 2i^d R - (1 + ^ 2/3)2 ) ct «p(^0
satisfies for all indices a, (3, j

\(Rd R r(qd q Y j e(q,R)\ < c aM {R)S~ j
To conclude the proof it remains to solve the equation for the difference <7\ =

-d RR - 2i^d R - (1 + fl2/3)2 j <Ji(R^,R) = e(R£*,R)

with zero Cauchy data at infinity. We claim that the solution <j\ satisfies

KRdnnqdgfa^R)] < c aM q-i{R)-\ j > 2

Note that this finishes the proof by defining a = a ap — o\. A change of vari-
able allows us to switch from the pair of operators (Rd R ,qd q ) to (Rd R ,£,d^) with
comparable bounds. We rewrite the above equation as a first order system for
{vi,v 2 ) = {(J 1 ,Rd R cj 1 ):

., 1 Vx \ j R^ 1 \ f Vl \ (

Cr

V2 J \ (l+R 2 /3) 2

Then we have

R- 1 - 2i£* j\v 2 )-\Re

d

-j-\ v f>-R-i\v\ 2 -R\v\\e

22

J. KRIEGER, W. SCHLAG, AND D. TATARU

which gives
and by Gronwall

^vl^-CiR-'M+Rlel)

POO / T}/ \ C

\v(R)\< J r (^J R'\e(R>)\dR>
Then for large j we have

(4.10) |e| < C i R~ 3 {R)~ i M < C i Rr°{R)~ 2 < q~ j (R}~ 2

To estimate derivatives of v we commute them with the operator. For derivatives
with respect to R we have

d R (Rd R ) ( M - ( U 5K , \ , ) (Rd R ) ( V ;

I

— 10J?.(l--R 2 /3)

V (l + fi 2 /3) 3

But the right-hand side is bounded by from the previous step and the hy-

pothesis on e, therefore as above ROrv is bounded by Rr' J .
We argue similarly for the £ derivatives. We have

. u 2

Z8 6 (Re)

The only difference is in the first term on the right, for which we write £5 =
and we use the decay property of v with j replaced by j + 1 :

\^V2\ < ^g-'- 1 < R-\~\ m(Re)\ < R-V

as desired. Finally, higher order derivatives are estimated by induction using the
above arguments at each step. □

Next, we describe the spectral measure of C. Due to the resonance at zero
energy, the spectral density becomes unbounded for small £. In what follows,
f-{R,0 :=f + (R,-0 = U(R,0-
Lemma 4.6. For all £ > there is a(£) ^ so that

<j>(R,Z)=a(Z)f + (R,Z) + a(Of-(R,0

with

i few ~ J 1 as £ ->

[ £ 2 a S £ — > oo

wifft symbol type behavior of all its derivatives. The density p(£) of the spectral
measure satisfies

' £-2 as £^0
£2 as £ — > 00
with symbol type behavior of all derivatives.

p(0

SLOW BLOW-UP SOLUTIONS

23

Proof. By inspection,

a{0 = w(f + U)J-U)) = Zt& w M->Q>U->®

By the preceding asymptotic analysis we can evaluate
with some small fixed e > to conclude that

£2 as £^0
as £ — > oo

with the corresponding upper bound on the derivatives. This yields the desired
upper bound on |a(£)|. To obtain the lower bound, we proceed as follows. First,
observe that

lm(/ + (i?,0/i(#,0) = -^

Second, it follows that

Im (f_(R, QW(f+(R, 0, 4>(R, 0)) =

so that

|Ty(/+(fl,0,W0)l>^ | l ^ ( ( ^

^From our asymptotic analysis, again at R — e£,~^ ,

W(R,0\ > / 1 as £^0

~U" 5 as £^oo

which leads to the claimed lower bound on |W(</>(-,£),
The Weyl solution

V>+ {R, £ + iO) = + m(f + iO)(j)(R, £)
satisfies tp + (-,£ + iO) = co (£)/+(■, £). Solving for c (£) yields

m ^ ] w(/ + (-,o,0(-,O) |iy(/ + (-,o,0(-,O)l 2

Since

/ + W) = -^,^(/ + (',9,9W)) + «WW/ + (',0>^6)

implies that

-2i£* = 0, /- (-, 0) = -2* Im [W(6(; £),/+(•, (-, 0, 0)]

we conclude that

(4.11) ,(0 - i ImmK + ffl) - -

The denominator was estimated above, leading to the desired bound on the spectral
density. □

24

J. KRIEGER, W. SCHLAG, AND D. TATARU

5. The transference identity

Returning to the radiation part e in ()3.5p . the idea is to expand it in terms of
the generalized Fourier basis <f>(R,£) from Proposition ^. 31 i.e., write

poo

i(r,R) = x (t)MR) + / x(t,OHR,Op(0^
Jo

and deduce a transport equation for the Fourier coefficients (xq(t),x(t,£)). The
main difficulty in doing this is caused by the operator R8r which is not diagonal in
the Fourier basis. We re-express this derivative in terms of the derivative 2£c>£. We
refer to this procedure, which involves a certain error operator, as the transference
identity since it allows us to transfer derivatives from R to £. We define the error
operator K, by

(5.1) Rd^u = -2£d 6 u + JCu

where / = Tf is the "distorted Fourier transform" from Proposition 14.31 and the
operator —2^8^ acts only on the continuous part of the spectrum. Apriori we have

AC : Cg°({&} U (0, oo)) -> C°°({&} U (0, oo))

Splitting functions on spec(£) into a discrete and continuous component we obtain
a matrix representation for /C,

K.

K*cd fccc

Using the expressions for the direct and inverse Fourier transform in Proposition ^. 31
we obtain

K-dd —

[R8r

MR), MR)) 2

1 L R

K-dcf =

a:

f(0Rd R <t>(R,0p(0<%, MR)) ,

' L ~R

ICcd(v) =

(r8r

MHR,v)) r2

' L R

<f

f(£)Rd R <KR,t)p{t)dt, <f>(R,r))) 2

' L R

■ (

(jT

2^f(OHR,Op(OdC,^(R,v)) T2

1 L R

ICdd = -L JCdcf = - I f(0K d (0p(0d^ lC cd (v) = K d ( V )

(5.2)

Integrating by parts with respect to R in the first two relations we obtain
1

2'

where

K d {r,) = {Rd R M4>{R,r } )
Integrating by parts with respect to £ in (|5.2[) yields

JCccfiv) = ( J KOlRdn - 2£dt]<l>(R,t)p(t)dt,<l>(R,Ti)

A + "m)^

SLOW BLOW-UP SOLUTIONS

25

where the scalar product is to be interpreted in the principal value sense with

/ecg°((o,oo)).

In this section, we study the boundedness properties of the operator JC. We
begin with a description of the function and of the kernel Kq(ji, £) of K, cc .

Theorem 5.1. a) The operator IC CC can be written as

(5.4) k; cc = -(^ + ^M)^- 7? )+/Co

where the operator K-o has a kernel Ko(r], £) of the forn^
(5-5) Ko (r,,0 = £ ^F(tv)

with a symmetric function F(£,rj) of class C 2 in (0, oo) x (0,oo) satisfying the
bounds

\d e F(£,ri)\ + \d v F(Z,r,)\ <
sup Wd*F(t,T,)\<

j+fc=2

where N is an arbitrary large integer.

b) The function is smooth and rapidly decaying at infinity.

Proof. We first establish the off-diagonal behavior of /C cc , and later return to the
issue of identifying the ^-measure that sits on the diagonal. We begin with (|5.3p
with / e C£°((0,oo)). The integral

p OO

<R) = / f(0[Rd R - 2Zd s ]<P(R,0p(0 d£
Jo

behaves like R at and is a Schwartz function at infinity. The second factor (j)(R, rf)
in (|5.3[) also decays like R at but at infinity it is only bounded with bounded
derivatives. Then the following integration by parts is justified:

vlCccf(v) = (u,£4>(R,vj) = (Cu,ct>{R^)) L2

+ ?/

-?? < 1

(t-

_1 (n

-?? > 1

i

C-

fry < 1

(e-

\-vY

-ni-

-r)i\

)" JV £-

1-77 > 1

ce +

-T) < 1

«-

br})'

" 2 (H

-77 > 1

Moreover,

/» CO />CO

£u=/ f(0[C,Rd R ]<l>(R,0p(0(%+ f(0(Rd R -2^d e )^(R,0p(0^
Jo Jo

I /"CO

with the commutator

2 The kernel below is interpreted in the principal value sense

26

J. KRIEGER, W. SCHLAG, AND D. TATARU

Thus,

Jo

Hence we obtain

V^ccfiv) ~ £cc(£/)(f7)

I ■ - double integral on the right-hand side is absolutely convergent, therefore we
can change the order of integration to obtain

(V - OK (V, = P (0(U(R)<P(R, £), <p(R, v)

mO(R9R-2^)4>(R,Op(0 dt
f(0U(R)4>(R,0p(0dC,cl>(R,v))

This leads to the representation in (15. 5ft when £ ^ 77 with
F(£,r)) = (u(R)4,(R,0,<t>(R,r))]

It remains to study its size and regularity. First, due to our pointwise bound from
the previous section,

(5.6)

sup|0(iU)

_R>0

\Rd R <f>(R,0

\d 6 <P(R,0

\d^(R,0

\d!<t>(R,0
\dl<p(R,0

we always have the estimates

<

< R V£ > 1

< min(i?^- 1 ,i? 3 )

< mm(RCKR 2 )

< mm{R 2 cKR 5 )
<mm(R 2 C\R 4 )

V£ > 1
VO < £ < 1

v<e > 1

V < £ < 1

(5.7)

\%F(Z,ri)

<rV x

< 4 2^2

< ^ 2 f] 2

V£ > 1, V > 1

> 1, 77 > 1

> 1, 77 > 1

They are only useful when £ and 77 are very close. To improve on them, we consider
two cases:

Case 1: 1 < £ + 77. To capture the cancelations when £ and 77 are separated we
resort to another integration by parts,
(5.8)

vFfotj) = (u(R)4>(R,0,jCct>(R,v)) = ([C,U(R)}<p(R,0,(b(R,v)) +£F(Z,v)

Hence, evaluating the commutator,

(5.9) (T)-tMt,v) =-((2U R d R + U RR )<t>(R, 0,HR,V)

Since U R (0) = it follows that (2U R d R + U RR )<f>(R,£) has the same behavior as
4>(R, £) in the first region. Then we can repeat the argument above to obtain

(V - 2 F(^v) = -([C,2U R d R + U RR ]<f>(R,0,<l>(R,ri))

SLOW BLOW-UP SOLUTIONS

27

The second commutator has the form, with V(R) := — 5(1 + i? 2 /3) 2 ,

[C, 2U R d R + U RR ] = 4[/ flfl £ - Wounds. ~ U RRRR - 2U R V R - 4tf. RH V r
Since V(R), U(R) are even, this leads to

(v - 2 Ht v) = ((c/ odd (i?)a fl + [/—(i?) + tu even (R))<j>(R, v)

where by JJ odd , respectively JJ even , we have generically denoted odd, respectively
even, nonsingular rational functions with good decay at infinity. Inductively, one
now verifies the identity

fc-l k

in , (v-0 2k mv) = ((E^'^f dR + J2^ en )^,0,H;V))
(5-10) j=0 1=0

(R)mf d (R)\ + \U e k T n (R)\ < {R)- 4 - 2k Vj,£
By means of the pointwise bounds on (f> and 8r4> from (|5.6p we infer from this that

( t) k -Hv)-*

(v - 2k

Combining this estimate with (|5.7|) yields, for arbitrary N,

£ (f + *7) _1 (l + l£* ~V*\)~ N provided f + Tj>l,

as claimed. For the derivatives of F we follow a similar procedure. If £ and ?y are
comparable, then from (|5 .T[) . \d v F(£, rj)\ < Otherwise we differentiate with

respect to 77 in (|5.10[) . This yields

fc-l fc

j=0 £=0

Using also the bound on F from above we obtain

fe-1

respectively

iw^)i< ^ ,, 2fc ^«i<e

and

(v-0 2

which again yield the desired bounds. Finally, we consider the second order deriva-
tives with respect to £ and r\. For £ and r\ close we again use the bound from (|5.7p .
Otherwise we differentiate twice in (|5.10[) and continue as before. We note that it is
important here that the decay of U£ dd and U^ en improves with k. This is because
the second order derivative bound at has a sizeable growth at infinity which has
to be canceled,

\%<KR,0)\ «i? 4

28

J. KRIEGER, W. SCHLAG, AND D. TATARU

Case 2: £, 77 <C 1. First, we note that -F(0, 0) = 0. This can be verified by direct
integration, and is heuristically justified by the fact that U = [C, ROr) — 2£. The
pointwise bound

sup \d e F^, v )\<i

t,v>o

follows by differentiating ()5.8|) and from the bound \d^<p(R,^)\ < R 2 , see (|5.6[) .
To bound the second order derivatives of F we recall the pointwise bounds, for

o<e<i,

\dU(R,0\

<

i? 2j R <

3 =0,1,2

If rj < £ < 1, then these bounds imply that

IVK,l)|< (R)- i R i dR+ i (R)- i R 3 C^dR+ ( (iZ)" 2 ^)-* dR

JO J£ 7 Jr; _ 2

<[i+iog(e/»?)]r'

The logarithm in the middle integral is an artefact and can be removed using the
oscillatory nature of dt<f)(R, £) in the regime R 2 £ > 1 as provided by Proposition ^. 51
and Lemma 14.61 Loosely speaking, this means integrating by parts using that

fy<t>(R,0 ~ R^dRe 1 ^ for R 2 £_ > 1 and small f. Thus, actually

\d iv m,v)\<r^

A similar computation yields, for £ < 1,

_ i

\d!F(Z,v)\< / (R)- 4 R i dR+ / ^Rr^R^dRZC*

JO J{~2

This bound is too weak when £ -C 77 < 1. In that case, we differentiate (|5.9|) to
obtain

(77 - OdlF(t 77) = 2%F(£, 77) + (dt<KR, 0, (2^ fl a fl + C/ flJi )0(i?, 77)

which in turn yields
(5.11)

-)2

(v-0%F(Z,r,) = J [2%F(C,C)
Using also the bound

\RdR^(R,C)\ <mm(R£-i,R 2 )
we can evaluate the inner product in (|5.1ip as follows:

0,{2U R d R + U RR )d ( <l>(;Q) d(

<)7

,Z),(2U R d R + U RR )d C <l>(;{)

ft 5 ri ' r°o

< / (R)- 6 R i R 2 dR+ / i (R)- 6 R 4 R(-z dR+ / i (R)- 6 R 2 ^- 1 RC" dR

JO JQ~^

<[i+iog(c/o]r'

cr 3

SLOW BLOW-UP SOLUTIONS

2!)

< 77*

The logarithm appearing in the middle integral can be removed as before exploiting
cancellations. Thus, (I5.11|) is controlled by

\(r,-OdlF(£,r,)\<\J\- h d{

Since 77 3> £ this yields

\%F(£,r,)\ < V -i

which concludes the analysis of the off-diagonal part of the kernel.

Next, we extract the 5 measure that sits on the diagonal of the kernel of JC CC
from the representation formula (|5.3p . see also (|5.4p . To do so, we can restrict £,77
to a compact subset of (0,oo). This is convenient, as we then have the following
asymptotics of <j>(R,£) for R£z ^> 1:

cf>(R,0=2Re [a{0e m ^]+O(R- 2 )

(Rd R - 2^)cf>(R,0 = -4 Re [& (Qe* 1 ** ] + 0(R- 2 )

where the O(-) terms depend on the choice of the compact subset. The R~ 2 terms
are integrable so they contribute a bounded kernel to the inner product in (|5.3[) .
The same applies to the contribution of a bounded R region. Using the above ex-
pansions, we conclude that the <5-measure contribution of the inner product in ([57
can only come from one of the following integrals:

(5.12) -4/ / f(Ox(R)^[^'(0<v)e iR(r2+r, ' 2) ]p(0^dR

Jo Jo

r°° r°° , 1 i.

(5.13) -4/ / f($ x (R)B*[Sa\$a(rj)e iR W^]p(QdedR

where x is a smooth cutoff function which equals near R = and 1 near R = 00.
In all of the above integrals we can argue as in the proof of the classical Fourier
inversion formula to change the order of integration. Integrating by parts in the first
integral (|5. 12[) reveals that it cannot contribute a ^-measure. On the other hand,
(|5.13|) contributes both a Hilbert transform type kernel as well as a <5-measure to
K. By inspection, the 6 contribution is

Re[£a'(£)a(r7)e iR «*- 7) *)]p(£) dR

)

= -87r^Re[ea / (Oo(0]p(0*tt-»?)

^ p(0 -

where we used that piO 1 = 47r£5| a (£)| 2 in the final step, see (|4.11[) . Combining

8(t - 77)

this with the <5-measure in (|5.3p yields (|5.4|) .
b) Arguing as in part (a) we have

KM - ^

For F we use the representation in (|5. 10|) with £ replaced by £d and replaced
by <fid- The conclusion easily follows from pointwise bounds on <fi(-,T)) and its
derivatives. □

30 J. KRIEGER, W. SCHLAG, AND D. TATARU

Next we consider the L 2 mapping properties for /C. We introduce the weighted
L 2 spaces L 2 ' a of functions on spec(£) with norm

/>oo

(5-14) := |/(^)| 2 + / \f(0\ 2 (® 2a p(0<%

Jo

Then we have

Proposition 5.2. a) The operators Kq, K, map

/C : L 2 p ' a -> L 2 p ' a+1 / 2 , /C : L 2 p a -> L 2 p ' a .
b) In addition, we have the commutator bound

: L 2 p a ^L 2 p a

with acting only on the continuous spectrum. Both statements hold for all

Proof. We commence with the /Co part, a) The first property is equivalent to
showing that the kernel

pHv)(v) a+1/2 K (v,mr a p-ko ■■ l 2 (r+)^l 2 (r+)

With the notation of the previous theorem, the kernel on the left-hand side is

We first separate the diagonal and off-diagonal behavior of Kq, considering several
cases.

Case 1: (£,7?) e Q := [0,4] x [0,4].

We cover the unit interval with dyadic subintervals Ij — [2 3 ' -1 , 2?' +1 ]. We cover
the diagonal with the union of squares

2

j=-oo

and divide the kernel Kq into

1qK = IahqKq + Iq\aKq

Case 1(a): Here we show that the diagonal part IadqKq of Kq maps L 2 to L 2 .
By orthogonality it suffices to restrict ourselves to a single square Ij x Ij . We recall
the Tl theorem for Calderon-Zygmund operators, see page 293 in [T^]: suppose the
kernel K(r], £) on M 2 defines an operator T : S — ► S' and has the following pointwise
properties with some 7 € (0, 1] and a constant Co:

(i) 1^(77,^)1 < Col^ - ^l" 1

(ii) \K( V ,0 - K(rf,S)\ <C \v~ VV\Z - vr 1 ^ for all |r? - rf\ < \£ - t?|/2

(in) \K(r,, - K(n, o\ < coie - mz - vr 1 - 1 for a u ^ - e\ <\c- W2

If in addition T has the restricted I? boundedness property, i.e., for all r >
and £0,% € R, ||T(oM>)|| 2 < C* r^ and ||T*(w^''°)ll2 < C r^ where w r '«°(£) =
— £o)A) w hh a fixed bump- function ui, then T and T* are L 2 (R) bounded
with an operator norm that only depends on Co-

Within the square Ij x Ij , Theorem 15.11 shows that the kernel of Kq satisfies
these properties with 7=1, and is thus bounded on L 2 .

SLOW BLOW-UP SOLUTIONS

31

Case 1(b): Consider now the off-diagonal part \q\aKq. In this region, by
Theorem 15.11

lA-ofaoi <(^r*

which is a Hilbert-Schmidt kernel on Q and thus L 2 bounded.

Case 2: (£,??) £ Q c . We cover the diagonal with the union of squares

B = |J li x I 3

3=1

and divide the kernel Kq into

1qcK = IbhqcKo + 1qc\ B Ko
Case 2a: Here we consider the estimate on B. As in case la) above, we use
Calderon-Zygmund theory. Evidently, |-ftTo(?7>£)l S l£ — ^l -1 on B by Theorem l5.ll
To check (ii) and (iii), we differentiate Kq, It will suffice to consider the case where
the derivative falls on F(^,rj). We distinguish two cases: if — 772 | < 1, then
|£ — t]\ < £2 which implies that

if, on the other hand, |£3 — 77^ | > 1, then

^1^4 V|£-£'|<K- 7/1/2
|^-77||^-»ji| |£"??| 2

which proves property (iii) on B with 7 = A and by symmetry also (ii). The
restricted L 2 property follows form the cancelation in the kernel and the previous
bounds on the kernel. Hence, Kq is L 2 bounded on B.

Case 2b: Finally, in the exterior region Q c \ B we have the bound, with arbi-
trarily large N,

\K (v,O\<a+O- N ^ + v)- N
which is L 2 bounded by Schur's lemma.

This concludes the proof of the first mapping property in part (a). The second
one follows in a straightforward manner since Kd is rapidly decaying at 00.

b) A direct computation shows that the kernel Kq"" 1 of the commutator [£c\ , Kq]
is given by

K° m (v,0 = (r,d v + Zd s )K (r,,Z) + Ko(ri,0 = 0-F com (H,r,)
interpreted in the principal value sense and with F com given by

pcom^ v) = t^-Fit, 7?) + + TfaWt, V )

By Theorem 15.11 this satisfies the same pointwise off-diagonal bounds as F. Near
the diagonal the bounds for F com and its derivatives are wors^l than those for F by
a factor of (1 = . Then the proof of the L 2 commutator bound for Kq is similar
to the argument in part (a).

The remaining part of the commutator [/C,£<9,e] involves

^The one derivative loss can be avoided by a more careful analysis, but this does not seem
necessary here.

:S2

J. KRIEGER, W. SCHLAG, AND D. TATARU

(i) The commutator of the diagonal part of /C cc with £e\. This is the multipli-
cation operator by

pb

which is bounded since p has symbol like behavior both at and at oo.

(ii) The operator S,d^lC c d which is given by the bounded rapidly decreasing func-
tion £0^(0-

(iii) The operator JC^d^ given by

!C dc Pd e f= K d {P)Pd^f{P)dP = - f(Od 6 m d (0)dP
Jo Jo

which is also bounded due to the properties of K~d- □

6. The final system of equations

We now rewrite the main equation (|3.5p in terms of the Fourier transform. With
T as in Proposition ^. 3[ and with j3 = AA _1 ,

which gives

f(d T + !3RdRj 2 = (d T + /3(-2^ + £)) T

= (d T - 2j3PdA 2 J^+20}c(d T - 20Pd^+ (3 2 (K. 2 + 2[£%,/C])JF
Recall that

/>oo

e(r, R) = x {T)<f> d + / x(t, P)<j>(R, P)p(0 dP
Jo

This leads to a transport type equation for the Fourier transform

X(t) = (x (t),x(t,O)

of £ by applying T to (|3.5[) . It is convenient to write it as a system for the two
components:

1 o (d T -2ppdA 2 + p ) x

(6 1)

= /3(J - 2K)(d T - 2(3Pd i )x - /3 2 (/C 2 - K. + 2[£% K])X

+ \~ 2 J r R(N 2k -i(R~ 1 J r ~ 1 X) + e 2k -i)

where it is understood that

N 2k „ 1 {R- 1 T- 1 X) + e 2k - 1 = (N 2k _ 1 (R- 1 T- 1 X) + e 2k _ 1 )(t(T),\- 1 (T)R)

Note that N 2k -\ and e 2k -\ are only defined on R < r, but for the Fourier transform
we need to extend them to all R - this will be described in the next section, but for
the moment just take an arbitrary compactly supported extension with the same
regularity.

We treat problem (|6.ip iteratively, as a small perturbation of the linear equation
governed by the operator on the left-hand side. For this we need to solve the

SLOW BLOW-UP SOLUTIONS

33

following uncoupled system consisting of an elliptic equation and a transport
equation:

d 2 T + Co

(d T - 2/3(r)^ s ) 2 + t] x(r, £) = 6(r, £),

(6.2)

We want to obtain solutions to (|6.1[) which decay as t — * oo. For the first equation
above this is achieved by using the standard fundamental solution Hq which has
kernel

-IWk-oi

This means that up to homogeneous solutions of the form e ^ d ^ T the unique
bounded solution to the elliptic equation is

1

Xd{r)

/ e -l^l 5 l T - CT l b d (a) da
Jo

The second equation is restricted to the range £ > 0. Hence it is hyperbolic in
nature, which means we can solve it backward in time, i.e., with zero Cauchy data
at t = oo. We denote by H the backward fundamental solution for the operator

[d T -2(3( T )^) +£
and by H(r,a) its kernel,

/oo
H{r,a)f{a) da

Combining the two components we obtain a fundamental solution for the system,

n = diag(H Q ,H)

Then we look for a solution X to (|6.1[) as a solution to the fixed point problem
X = h((3(I - 2K)(d T - 2j3£d^X - (3 2 {JC 2 - K + 2[^,K])X

(6.3)

\- 2 TR{N 2k -i{R' 1 T- 1 X) + e 2k -i)

Remark 6.1. One can also freely add Ce~^ d ' ,T to the first component x d of X.
Thus the fixed point argument yields in effect a one parameter family of solutions
X depending on the parameter C.

The mapping properties of H are described in the following result, which was
proven in [TP] ,

Proposition 6.2. For any a > there exists some (large) constant C = C(a) so
that the operator H(t, a) satisfies the bounds

( 'er\ c

(6.4)
(6.5)

uniformly in a > r

\\H{T,0)\\ , .:

(d T -2P(r)^)H(T,o)

<

34

J. KRIEGER, W. SCHLAG, AND D. TATARU

This leads us to introduce the spaces L°°' N L 2 p ' a with norm
\\f\\ LX , NLr :=supr N \\f(T)\\ L2 , a

P T>1 P

Then the above proposition immediately allows us to draw the following conclusions:
Corollary 6.3. Given a > 0, let N be large enough. Then

\\Hb\\ L ~,N-* Lr +^ + I (dr - 2P{T)Zd^Hb\ LaatN _^ a < C Q 1||6||

with a constant Co that depends on a but does not depend on N .

The small factor N 1 is crucial here for our argument to work. For Hq we have a
stronger straightforward counterpart of the above result:

Lemma 6.4. The operator Ho satisfies the bounds

||-ffoMlz,°°.- N + II^t-HoMIl 00 '™ — GivllMU 00 '^

We note that while the constant on the right cannot be small, we no longer lose
powers of t compared to the bounds for H . Hence for fixed ./V we can choose To
depending on N so that we gain the smallness in the range r > r . Combining the
last two results we obtain

Proposition 6.5. Given a > 0, let N be large enough. Then there exists tq
depending on N so that for t > tq we have

\\nb\\ L ^ N _ 2L2 , a+l/2 + (d T -2f3(T)^dAnb <Coh\b\\ LOO , NL 2,«

with a constant Co that depends on a but does not depend on N .

On the other hand, the nonlinear operator N2k-i from (|6.ip has the mapping
properties stated in Proposition 16.71 below. We first relate the spaces L 2 p ,a to the
Sobolev spaces in R 3 .

Lemma 6.6. Let a > 0. Then

||x|| L 2, Q x \\Rr 1 T~ 1 x\\ H z<*(m?)

Proof. For integer k we have

k k

\\x\\ L *.* x w&f-^h* ~ E wr-'&f-^Wl^)

j=0 j=0

k

3=0

But

R^CR = -d 2 R - ^d R - 5W\R)

where the first two terms can be recognized as the radial part of the three-dimensional
Laplacian. Hence we get

k

NIl;-* x £ II (- A - ^w\R)yR- 1 T- 1 x\\ L 2 m

3=0

SLOW BLOW-UP SOLUTIONS

35

Since W is bounded together with all its derivatives, the conclusion of the lemma
follows for integer a.

For noninteger a we use interpolation. First we consider the map

x i — * R~ x T~ x x

and obtain the bound

\\R~ 1 f~ 1x \\h^(«?) ^ IMIz,2.°

To obtain the reverse bound we use the map

u i— ► J-RS(u)

where S(u) stands for the spherical average of a function u in R 3 . □
Proposition 6.7. Assume that N is large enough and g < a < j. Then the map

x -» X- 2 TR{N 2 k-i{R~ X J r ~ X x))
is locally Lipschitz from L°°' N - 2 L 2 p ' a+1/2 to L°°' N L 2 ' a .
Proof. Using the lemma, it remains to prove that the map

v ^ \~ 2 N 2k -i(v)
is locally Lipschitz from L°°' N - 2 H 2a+x to L°° N H 2a . We recall that

N 2 k-i(v) = 5( u tk-i -Ua)v + lOui^-i v 2 + 10u^_ 1 v 3 + 5u 2k -i v 4 + v 5 ,

see the comment following (|3.1[) , The time decay is trivially obtained for all but
the first term, for which we need an additional step, where we pull out a factor of
u 2k -i — uq. Using the regularity of u 2k -i given in Theorcm l2.ll we obtain

A*

u 2k -i-u Q e ^2 IS 2 (R,Q k -i), A~ 2 ("Uq - u 2 k-i) e t- 2 S 2 (R- 2 , Qfc-i)

This indicates that two units of decay in r are gained. On the level of Sobolev
spaces we argue as follows: since we are working with inhomogeneous Sobolev
spaces, we can localize the above estimate to unit cubes, as the £ 2 summability for
N 2 k-i(v) with respect to unit cubes is inherited from any of the v factors. But
in any unit cube Q the coefficients u 2 k-i have at most (1 — a)~^~~ singularities
(where a = j ~ -p) so that we can bound them in Sobolev spaces

A 2 A 2

\\u2k-i\\m+2«(Q) < dist (Q Q y \\U2k-i - uo\\m+^(Q) < ^dist(Q,0)
where we used that a < K. Then it suffices to establish the quintilincar estimate

jj2a+l _ jj2a+l _ jj2a+l jj2a + l jj2a + l ^ jj2a

which in three space dimensions holds for a > i (a standard application of the
fractional Leibnitz rule and Sobolev imbedding, see [22], page 105). □

3(5

J. KRIEGER, W. SCHLAG, AND D. TATARU

7. Conclusion

We now prove Theorem 11.11 We first construct a blow-up solution inside the
cone as follows. We begin with the approximate solution v,2k~i and the error e2k-i
given by Theorem l2.1l inside the cone. We extend them outside the cone to functions
having the same size and regularity, supported in r < 2t. Then the relation

2

&2k-i = (~dt + d 2 + ~d r )u2k-\ + ulk-i
r

is valid only inside the cone.

Using Propositions 16.51 16.71 we iteratively find a solution

X e L°°> N - 2 L 2 p > a+1 l 2

for the equation (|6.3p for t < to sufficiently small and

1 v

Then we set

v = R-^-'X G L°°< N - 2 H 2a+1

and

u = U 2 k-1 + v
Given the derivation of (|6.3p . the function u solves

2 2
{-d 2 + d 2 r + -d r )u + u b = {-d 2 t +d 2 + -d r )u 2k -i + ~ e 2k -i
r r

which implies that the function u solves the nonlinear wave equation

(-d? + a 2 + -d r )u + u 5 = o

r

inside the cone.

The second part of the argument is to extend the above solution u to the exterior
of the cone K = {0 < t < t , 0< |x| <£}so that the blow-up occurs only at the
tip of the cone. For this we first observe that the above function u is close to u
inside the cone and close to outside, namely

lim / \V(u(t) - u (t))\ 2 + \u(t) -u (t)\ 6 dx =

and

limy \Vu(t)\ 2 + \u{t)\ 6 dx =

Hence given 6 > we can choose to so that the two quantities above are less than <5 6 .
We let w be the solution to the equation

{-d 2 + d 2 + -d r )w + w 5 =
r

with initial data

w(t ) = u(t ), w t (t ) = u t (t Q )

Due to the finite speed of propagation we conclude that w = u inside the cone. To
conclude the proof of the theorem we will show that w cannot blow up outside the

SLOW BLOW-UP SOLUTIONS

37

cone before time 0. For this it suffices to prove that the energy outside the cone
stays small,

(7.1) / \Vw(t)\ 2 + \w(t)\ 6 dx < S

see Q2], Q1-

This is proved using energy conservation. The energy £(w(t)) is conserved in
time. At time to we have

£(w(t )) = £{u(t ))=S(W) + O(S)

hence at time t we must have a similar relation. But the energy inside the cone is
already close to this, so we obtain

On the other hand, we have the Sobolev inequality

< S

M 6 dx <

AT

with a universal constant (independent of t) . Combining the two inequalities above
we see that for each t there are two possibilities. Either we have

or

f \{ut + \Vu\Ydx<5
I \{ul + \Vu\fdx>\

'at

The first alternative holds at t = to. Then a continuity argument shows that it
must hold at all t, since the above integral is a continuous function of t for as long
as it stays small.

References

[1] Aubin, T. Nonlinear analysis on manifolds. Monge-Ampere equations. Grundlchren der
Mathcmatischen Wissenschaften, 252. Springer- Verlag, New York, 1982.

[2] Bizori, P., Chmaj, T., Tabor, Z. On blowup for serailinear wave equations with a focusing
nonlinearity. Nonlinearity 17 (2004), no. 6, 2187-2201.

[3] Dunford, N., Schwartz, J. Linear operators. Part II. Wiley Classics Library. John Wiley &
Sons, Inc., New York, 1988.

[4] Gesztesy, P., Zinchenko, M. On spectral theory for Schrodinger operators with strongly sin-
gular potentials, Math. Nachr. 279 (2006), 1041-1082.

[5] Grillakis, M. Regularity and asymptotic behaviour of the wave equation with a critical non-
linearity. Ann. of Math. (2) 132 (1990), no. 3, 485-509.

[6] Jorgens, K. Das Anfangswertproblem im Grossen fur eine Klasse nichtlinearer Wellen-
gleichungen. Math. Z. 77 (1961), 295-308.

[7] Karageorgis, P., Strauss, W. Instability of steady states for nonlinear wave and heat equations,
preprint 2006.

[8] Kenig, C, Merle, F. Global well-posedness, scattering and blow-up for the energy critical

focusing non-linear wave equation, preprint 2006.
[9] Krieger, J., Schlag, W. On the focusing critical semi-linear wave equation, preprint 2005, to

appear in Amer. Journal Math.
[10] Krieger, J., Schlag, W., Tataru, D. Renormalization and blow up for charge one equivariant

critical wave maps., preprint 2006.

38

J. KRIEGER, W. SCHLAG, AND D. TATARU

[11] Levine, H. Instability and nonexistence of global solutions to nonlinear wave equations of the

form Putt = -Au + T(u). Trans. Amer. Math. Soc. 192 (1974), 1-21.
[12] Marchcnko, V. A. Sturm- Liouville Operators and applications. Operator Theory: Advances

and Applications, Vol. 22, Birkhauser 1986.
[13] Merle, F., Zaag, H. Determination of the blow-up rate for the semilinear wave equation.

Amer. J. Math. 125 (2003), no. 5, 1147-1164.
[14] Merle, F., Zaag, H. On growth rate near the blowup surface for semilinear wave equations.

Int. Math. Res. Not. (2005), no. 19, 1127-1155.
[15] Merle, F., Zaag, H. Determination of the blow-up rate for a critical semilinear wave equation.

Math. Ann. 331 (2005), no. 2, 395-416.
[16] Schlag, W. Spectral theory and nonlinear partial differential equations: a survey. Discrete

Contin. Dyn. Syst. 15 (2006), no. 3, 703-723.
[17] Shatah, J., Struwe, M. Geometric wave equations. Courant Lecture Notes in Mathematics,

2. New York University, Courant Institute of Mathematical Sciences, New York; American

Mathematical Society, Providence, RI, 1998.
[18] Sogge. C. D. Lectures on nonlinear wave equations. Monographs in Analysis, II. International

Press, Boston, MA, 1995.
[19] Stein, E. Harmonic Analysis, Princeton, 1994.

[20] Strauss, W. Nonlinear wave equations. CBMS Regional Conference Series in Mathematics,
73. AMS, Providence, RI, 1989.

[21] Struwe, M. Globally regular solutions to the u 5 Klein-Gordon equation. Ann. Scuola Norm.
Sup. Pisa CI. Sci. (4) 15 (1988), no. 3, 495-513 (1989).

[22] Taylor, M. Tools for PDE. Mathematical Surveys and Monographs, 81. American Mathe-
matical Society, Providence, RI, 2000.

Harvard University, Dept. of Mathematics, Science Center, 1 Oxford Street, Cam-
bridge, MA 02138, U.S.A.

E-mail address: jkriegerOmath.harvard.edu

Department of Mathematics, The University of Chicago, 5734 South University Av-
enue, Chicago, IL 60637, U.S.A.

E-mail address: schlagSmath.uchicago.edu

Department of Mathematics, The University of California at Berkeley, Evans Hall,
Berkeley, CA 94720, U.S.A.

E-mail address: tataruOmath.berkeley.edu